Self Force Orbit-Integrated gravitational waveforms for E(I)MRIs 
(Talk given at the 9*^^ LISA Symposium, 21-25 May 2012, BnF-Paris) 
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We present the first orbit-integrated self force effects for an IMRI or EMRI source, specifically 
the effects of its conservative piece on the orbit and on the waveform. We consider the quasi- 
circular motion of a particle in the spacetime of a Schwarzschild black hole, find the orbit and the 
corresponding gravitational waveform, and discuss the importance of the conservative piece of the 
self force in detection and parameter estimation. We also show the effect of the conservative piece of 
the self force on gauge invariant quantities, specifically u' as a function of the angular frequency Q. 
For long templates the inclusion of the conservative piece is crucial for gravitational-wave astronomy, 
yet may be ignored for short templates with little effect on detection rate. 

An important source of gravitational waves for space borne detectors are E(I)MRIs, extreme (intermediate) mass 
ratio inspirals. Such sources will allow us to test directly the Kerr hypothesis, and allow us to map the spacetime 
surrounding black holes. Moreover, the detection of E(I)MRIs will allow us to determine the mechanisms that shape 
stellar dynamics in galactic nuclei with unprecedented precision ([l| and references cited therein). 

The orbits of E(I)MRIs are typically highly relativistic, and exhibit exciting phenomena, e.g. extreme periastron 
and orbital plane precessions. Because the orbital evolution time scale ("radiation reaction time scale") is much 
longer than the orbital period(s), over short time scales the orbit is approximately geodesic, yet on long time scales 
it deviates strongly from geodesic motion of the background. Instead, the smaller objects moves along a geodesic 
of a perturbed spacetime. Equivalently, one may construe the orbit as an accelerated, non-geodesic motion in the 
spacetime of the unperturbed central object, where the acceleration is caused by the self force of the smaller object 
(0). We undertake here the latter approach. 

We consider here a point particle of mass /i that moves in a quasi circular orbit (i.e., an orbit that would be circular 
but for the radiation reaction effects) in the space-time of a Schwarzschild black hole of mass M , under the assumption 
that fi ^ M. Specifically, we take here fi = 10~^M, and the orbit starts at rg = lOM down to close to the ISCO 
at r = 6 M. The orbital evolution is driven by the particle's self force. In practice we use the self force found for 
circular orbits in the Lorenz gauge and interpolate (using the least number of necessary terms) the tabulated values 
to obtain smooth functions for the self force to the accuracy given in [2]- Specifically, we expand the components of 
the self force as follows: 
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where PN5.5 stands for the standard -^-post-Newtonian expression (converting Eq. (3.1) in \7] from luminosity to /*). 
Fitting the free parameters, we find the values appearing in Table 1. Finally, the last remaining non-zero component 



of the self force (w*") is found from the condition that the self force must satisfy /, 
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0, where the 4-velocity is 



taken for the circular geodesic orbit for which the self force was calculated. 

We evolve the orbit using three different codes. First, we use the energy balance approach to evolve the orbit. This 
approach ignores conservative effects on the orbital evolution and correspondingly on the waveforms, and can be done 
without the self force: fluxes of otherwise conserved qualities to infinity and down the event horizon are used to update 
the particle's constants of motion. As expected, we show in [1| that the orbit and the resulting waveforms are identical 
in the energy balance approach and in our self force calculations, when the conservative piece of the self force is turned 
off (i.e., when we take bf = 0.) We used two independent codes for the computation of the orbital evolution with the 



2 



TABLE L The fit parameters for tfie self force. These parameters reproduce the accuracy of 01 to all significant figures for 
all data points. Our results for a^gj^ are very crude predictions for the corresponding PN parameters, as our fit ignores all 
higher-order terms. 
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FIG. 2: The osculating-code variables a and /3 as functions of the time t (upper two panels), and the effective eccentricity e of 
the orbit (lower panel) as a function of t in the osculating case. 



3 



self force. First, we used the method of osculating geodesies (specifically eqs. (43)-(47) in @). We also integrate the 
equations of motion directly, i.e., solve directly for u'^ from the equation of motion m/^V^m" = Ai"VsF ("Newton's 
second law," with covariant differentiation compatible with the background metric) and integrate its solution to find 
x^. Both codes are numerically stable and convergent. Specifically, the osculating code converges with 5*^^ order, and 
the direct code converges with 4"^ order (Fig. [T]). 

The integration using the method of osculating geodesies cannot keep the value of the eccentricity as precisely zero. 
As both variables a and /3 (see 6] for definitions) are dynamical, the eccentricity e must evolve along the orbit too. 
This behavior is shown in Fig. [21 Interestingly, the inclusion of the conservative piece of the self force amplifies the 
resulting effective eccentricity. 
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FIG. 3: The Orbit. The 4-position for the three orbital evolution codes: energy balance (dotted), direct evolution (dashed), 
and osculating code (solid). 
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FIG. 4: The Orbit. The 4-velocity for the three orbital evolution codes: energy balance (dotted), direct evolution (dashed), 
and osculating code (solid). 

The orbit is displayed in Figs. [3] and |4] and Figs. [5]|6] for the three codes. Notably, the two independent self force 
codes reproduce the orbit to high level of agreement, with a difference much smaller than the difference between either 
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FIG. 5: The Orbit. The shape of the orbit for the three orbital evolution codes: energy balance (dotted), direct evolution 
(dashed), and osculating code (solid). 
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FIG. 6: The Orbit. A small portion of the orbit is magnified to show detail. 

and the orbit generated in the energy balance approach. This difference is attributed to the effect of the conservative 
piece of the self force. 
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The orbit, of course, is a gauge dependent quantity. Indeed, the position vector changes trivially under gauge 
transformations, x'^ — >■ + . We can, however, create gauge invariant quantities with a specific gauge choice (in 
our case, the Lorenz gauge), and then those quantities - by virtue of their gauge invariance - are guaranteed to remain 
unchanged in any other gauge. Two independent gauge invariant quantities are u* ( "gravitational redshift" ) and the 
angular frequency (3). In Fig. [7] we plot as a function of with and without the conservative piece of the self 
force. Notably, to the accuracy of our numerical computation the two curves overlap. That is, wc find that u* as a 
function of fl is insensitive to the conservative piece of the self force. This conclusion implies that when an actual 
data stream is used and this gauge invariant figure is plotted, one may use a simplified radiation-reaction scheme, 
that does not include the conservative effects in its analysis. 

There is, however, an aspect of the gauge invariant figure that is sensitive to the conservative effects, specifically 
the speed with which the data point moves along the curve. The way the conservative effects are manifested in the 
gauge invariant plot is not is the shape of the curve, but in the time it takes the signal to move along it. 
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FIG. 7: The gauge invariant figure of u* as a function of the angular frequency Q. The curve without (dashed, o) and with 
the conservative effects (solid, □) are shown, together with equal-t spacing marks starting at t = 1,500M in increments of 
1,000M. The two curves are indistinguishable to the numerical accuracy of our computation (using the osculating code in 
both cases). 

After we obtain the orbit we used it in order to evolve the linearized Einstein equations, using a code for the 
sourced Teukolsky equation with hyperboloidal slicing (Q). We present the waveforms in Fig. [Sj We show the 
three waveforms, for the energy balance approximation — which is equivalent to the self force waveform with the 
conservative piece turned off — and the two independent waveforms obtained from the direct approach and from the 
osculating geodesies approach that use the full self force, including its conservative piece. All three waveforms start 
in phase at early times. The two full self-force waveforms are dephased from the energy balance waveforms by much 
more than from each other. Specifically, the total commutative dephasing of either waveform from the energy balance 
waveform over the entire orbit from rg = 10 M down to near the ISCO, is A0 = 8.4 ± 0.4 rad. This dephasing 
corresponds to about 4/3 of a cycle, compared with the 54.3 cycles the particle makes in the energy balance case. 

We next study the importance of the conservative piece of the self force for the waveforms. Specifically, we take a 
window of duration L from the end of the chirp part of the energy balance waveform, and find its overlap integral 
with the full self force waveform, which in this case models the actual data stream. As expected, the longer L, the 
smaller the overlap integral. When we take L = ithrcshoid — 816.6 M we find that the overlap integral is reduced to 
0.96, which corresponds to a loss of 10% in the event rate. Therefore, for L < ^threshold ignoring the conservative 
piece of the self force does not lead to a significant loss of accuracy in detection of E(I)MRI events. However, longer 
stretches of data, L > ^threshold leads to a significant decrease in event rates, so that the inclusion of the conservative 
piece of the self force is important. Further detail of this work appears in [J]. 

The authors are indebted to Gaurav Khanna for discussions and for use of his numerical code. This work was 
supported by a NASA EPSCoR RID grant and by NSF grants PHY~0757344 and DUE-0941327. LMB is grateful to 
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FIG. 8: The real part of the waveform -04 when the conservative piece of the self force is included (solid for the osculating orbit 
method, and dashed for the direct integration method) , and using the energy balance approach (dotted) . The three waveforms 
are in phase at the beginning of the waveforms, but the former two dephase with time from the latter (and to a much smaller 
extent from each other). 



Alessandro Spallicci for hospitality. 
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